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We investigate the collective excitations of an atomic Bose- 
Einstein condensate in the self-binding regime produced by 
electromagnetically induced "gravity" (1/r attraction). An- 
alytical expressions for the frequencies of the monopole and 
quadrupole modes are obtained at zero temperature, using 
the sum-rule approach, and compared with the exact results 
available in the Thomas-Fermi limit. The low-energy dynam- 
ics of such condensates is shown to be dominated by the effec- 
tive "plasma" frequency. An analog of the Jeans gravitational 
instability is analyzed. 

PACS number(s): 03.75.Fi, 34.20. Cf, 34.80. Qb, 04.40.-b 

Introduction — Recently, a new kind of atomic Bose- 
Einstein condensates has been proposed H. It has 
been shown that particular configurations of intense, off- 
resonant electromagnetic fields mimic the gravitational 
attraction between atoms located well within the wave- 
length of these fields, giving rise to scenarios very dif- 
ferent from the ones characterizing usual Bose-Einstein 
condensates The interatomic potential is of the 

form —u/r, where the gravitational coupling u — Gm 2 is 
replaced by 
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a 2 q 2 I . 
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Here the total laser intensity and wavenumber are de- 
noted by / and q, respectively, and a is the atomic polar- 
izability at the frequency cq Q| . For a sufficiently strong 
induced "gravity" the BEC becomes self-bound, i.e. sta- 
ble in the absence of an external confining potential j| . 
Such a self-bound system may become a laboratory ana- 
log of a non-relativistic Bose star B, a self-gravitating 
system that balances gravity with the "quantum pres- 
sure" corresponding to the zero-point kinetic energy fluc- 
tuation. In the self-bound condensate considered in Ref. 
|jj the induced "gravity" can be also balanced by the 
short-range interatomic repulsion. 

It is known that the long-wavelength collective exci- 
tations of a spatially-uniform system of both fcrmions 
(e.g., the electron gas) and bosons jfj], interacting via a 
repulsive potential of the form u/r, are characterized by 
a spectral gap, fixed by the "plasma" frequency 



lu 2 = Aitun/m 



(2) 



where n is the particle density. In this paper we show 
that also in the case of a self-bound atomic BEC, inter- 
acting via an attractive —u/r potential, the collective 
oscillations are characterized by an effective "plasma" 



frequency io v of the form (Q), with a spatially- averaged 
density n and the positive coupling constant u [Eq. (yj)]. 
Its correspondence with the fundamental oscillation fre- 
quency of white dwarfs and neutron stars |t]] is pointed 
out and the analog of the Jeans gravitational instability 
is explicitly analyzed. 

Generalized Gross- Pitaevskii equation — The Gross- 
Pitaevskii equation ||] for the order parameter ^(r, t) 
can be obtained starting from 
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~dt ~ <54" 
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Here the total mean- field energy functional -fftot = -ffkin + 
H\ w + fZgrav + H s consists of: (a) the kinetic energy, 
H kin = j dr(h 2 /2m)\V^\ 2 ; (b) the harmonic-trap en- 
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ergy, H ho = JV ho ndr, where n = |*(r, t)\ 2 
atomic density and Vho = mui 2 r 2 /2 is an isotropic har- 
monic external potential; (c) the mean-field energy due 
to the electromagnetically induced "gravity", H grav — 
(1/2) J nQdr, where $ is the "gravitational potential" 
solution of the Poisson equation 



V~$ = Airun 



(4) 



with the boundary condition $ = for \r\ — * oo; (d) the 
mean-field contribution due to the short-range scattering, 
H s = (g/2) J n 2 dr, where g = Airah 2 /to, a being the s- 
wave scattering length. 

By expressing the complex BEC order parameter 
through the condensate density n and phase 4> as ^ = 
y/nexp(i(j>), we can write Eq. (||) as a coupled set of 
collisionless hydrodynamic equations: 



dn 
~dt 



V(nv) = 0, 



(5) 
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+ gn + V ho + $ = . (6) 



The last equation establishes the irrotational property 
of the superfluid velocity field defined by v = hV<fi/m. 
Similar equations have been derived in the case of har- 
monically trapped Bose-Einstein condensates interacting 
via zero-range forces || . 

Jeans-like instability — In the absence of external po- 
tentials, Eqs. (||) and (^|) admit a stationary solution 
with constant density n. Such a solution is, however, 
dynamically unstable, as can be seen directly from the 
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Bogoliubov dispersion relation JlTJ] corresponding to (||) 
and (i): 



h 2 k 4 
4m 2 



„2, 2 



(7) 



In (0) c s = y/gn/m is the Bogoliubov speed of sound 
and the plasma frequency u p is given by (|^). The oc- 
currence of imaginary frequencies LOk for small values of 
k reveals the existence of an instability, which resem- 
bles the Jeans gravitational instability . Assuming 
Lu p -C mc 2 , the minimal wavelength ensuring a real value 
for the frequency (0) is given by 



Aj = 2'Kc s /ujp — 2ir(h 2 a/mu) 1 / 2 



(8) 



In a "self-gravitating" BEC, both uj p and c s scale as n 1 / 2 , 
so that the critical wavelength Aj is independent of the 
density, as explicitly revealed by (^). 

Besides these unstable solutions, Eqs. (j|) and (^) ad- 
mit a stable self-bound solution for any finite number N of 
atoms. In this case long-wavelength fluctuations, which 
can trigger the Jeans instability, are excluded by the fi- 
nite size of the system. Remarkably, the ground-state 
density of a self-bound Bose cloud, in the limit in which 
the kinetic energy can be neglected (Thomas-Fermi limit 
for Bose gases), has the analytic profile Q 



n (r) 



N sin(7rr/i? ) 
4R 2 r 



6(R - r) = 0, 



(9) 



where the radius Rq — Xj/2 is fixed by the Jeans wave- 
length and 8 is the Heaviside function. 

Sum rule approach — The quadrupole and monopole 
excitation frequencies at zero temperature can be esti- 
mated using a sum-rule approach |j],[l2| . The calculation 
of the oscillation frequencies ui is based on the ratio 



h 2 u 2 = m^/m\ 



(10) 



between the third moment (cubic-energy weighted) and 
the first moment (energy- weighted) of the dynamic struc- 
ture factor: m p = J2„( huj „o) p |< | F \ n >| 2 (p = 1,3). 
Here hu> n Q is the excitation energy of the state |n > 
and F is an excitation operator chosen similarly to |)J: 

F = YliLi r i or F — J2?=i r i Y 2m, for the monopole 
or quadrupole operators, respectively. The moments mi 
and 77x3 can be reduced to the form of commutators in- 
volving the Hamiltonian of the system and evaluated in 
terms of ground-state expectation values. One then ob- 
tains the following results: 
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6H ho + 6H S 



Nm(r 2 ) 
4^ + 4^0-4/5 5, 
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(11) 
(12) 



where, in deriving the second equation ([11]), we have used 
the virial identity 
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The frequencies u) m and ujq , given by ([11]) and ( |12|) , arc 
always real for positive values of the scattering length. 
If the expectation values of Eqs. ( |j"l| ) and ([jj]) are eval- 
uated using the exact ground state, they provide a rig- 
orous upper bound to the lowest excitation frequencies 
for a given multipole. This approach has been shown to 
provide an excellent description of the crossover between 
the non-interacting and the Thomas- Fermi regimes in the 
case of condensates interacting via zero-range forces || . 

A useful analytical estimate for the expectation val- 
ues of Eqs. (Jll]) and ([l2]) can be obtained by replacing 
the exact ground state with the gaussian wavefunction 
^ = iVa exp (— r 2 /2w 2 )) /tt^w^ . Here w is a variational 
parameter that can be calculated by minimizing the to- 
tal energy H to t and is related to the expectation value 
(r 2 ) = (3/2) w 2 . The various contributions to the en- 

H kin = 3Nh 2 /Amw 2 , H ho 



I \l2i:mw i and H t 
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ergy H tot then take the form: 
= Nmuj 2 w 2 /A, H s = ah 2 N 2 

— {1/\/2tt)uN 2 /w. The gaussian ansatz allows us to de- 
scribe analytically the various regimes in the "phase di- 
agram" of the problem (!]]. It is noteworthy that the 
monopole and quadrupole frequencies obtained by the 
gaussian ansatz coincide with the ones derivable from a 
time-dependent variational calculation in a gaussian ba- 
sis @. 

Low-energy excitations in a self-bound BEC — In the 
following wc calculate the monopole and quadrupole ex- 
citation frequencies for a self-bound BEC from Eqs. ( |Tl] ) 
and (]l^), using, for the ground state, either the exact 
solution of the Gross-Pitaevskii equation or the gaussian 
approximation. In these calculations we are interested in 
the effects of s-wave scattering, measured by the dimen- 
sionless parameter 



s = muN 2 a/h 2 



(14) 



Upon neglecting the role of the external trap, the 
width of the ground state, within the gaussian ansatz, is 

given by w = (3/2)V?i72(l+ Vl+P 7r)s J A, where 

A = H 2 /muN is the effective "gravitational" radius which 
fixes the size of the condensate in the absence of s-wave 
interaction (a = 0). This expression for w allows for the 
analytical evaluation of the monopole and quadrupole 
frequencies ( |ll|) and (|lj) that are plotted in Fig 1 in units 
of the " gravitational" frequency mu 2 N 2 /h 3 . Comparison 
with the results obtained using the exact ground state 
shows that the gaussian ansatz provides indeed a good 
approximation, for both positive and negative values of 
s, provided one is far from the critical value s c ~ — 1, 
where the condensate exhibits collapse 

The order of magnitude of lum and ujq is given by the 
"plasma" frequency ui p , irrespective of the scattering pa- 
rameter s. This can be simply seen by using the virial 
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identity (|l|). Assuming s > and setting H\ 1Q = 
(absence of external trapping), we derive the inequalities 
< 2i7kin < —Hgnw and, hence, from ( O ) and ([TJ]), we 
obtain 1 < u) 2 m /uj 2 < 2 and (4/5) < uJq/u> 2 < (14/5). 
Here ui 2 = —H glav /Nm(r 2 } turns out to be always pro- 
portional to the plasma frequency oj 2 = Airun/m. Sim- 
ilarly, assuming s c < s < 0, we obtain < u> 2 4 /ui 2 < 1 
and (14/5) < oJq/oj 2 < (24/5). The proximity of the 
monopole frequency lum to the "plasma" frequency to p 
breaks down upon approaching the critical value s c ~ — 1 , 



where 



J M 



(whereas ujq — > (24/5) or). 
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FIG. 1. Monopole and quadrupole frequencies of a 
self-bound condensate (in mu 2 N 2 /Ti 3 units) versus the scat- 
tering parameter s = muN 2 a/h 2 , for both positive and neg- 
ative values of a. Circles: monopole frequencies; squares: 
quadrupole frequencies. The solid and long-dashed lines cor- 
respond to the gaussian approximation. Inset: the corre- 
sponding Log-Log plot for a > 0. 

The limits of large and small values of the scatter- 
ing parameter s (denoted as the TF-G and G self-bound 
regimes, respectively [Q) will now be discussed sepa- 
rately: 

(a) TF-G-regime — In the TF-G-regime (self-bound 
Thomas-Fermi "gravitational" regime) [|J "gravity" is 
balanced by the s-wave interaction (large- s limit), while 
both the kinetic energy fljdn (oc Ti 2 ), and the external 
trap potential are negligible. In this regime the monopole 
and quadrupole frequencies in (O ) and ( |l2| ) are given, 
using the ground state profile (Q)7 by ujm = 0.6225 oj p , 
luq = 0.3928 Lu p , where 



lu p = 0.5642 m*u*Nzh~ 



(15) 



is the plasma frequency evaluated at the central density. 
The ratio ujm/^q turns out to be equal to \/5/2. Re- 
markably, this is the same result holding for a harmon- 
ically trapped BEC in the Thomas-Fermi limit in the 
absence of the 1/r interaction 0. In the TF-G regime, 



the chemical potential [i = (H] i i n +Hi 10 +2H glav +2H s ) /N 



has the value — /i 



>> hu> r , 



It is interesting to note that in the Thomas-Fermi limit 
the sum-rule expression (11) for the monopole frequency 



reduces to 



= 2\H giav \/Nm(r 2 ^ 



and coincides with 



the value derived in politropic stars using classical hy- 
drodynamics j7|. 

In the TF-G regime the static linear response for 
the monopole operator F — Yli=x r 1 can be calcu- 
lated analytically. This allows for an alternative way 
to estimate the monopole frequency through the ratio 
h 2 uj 2 = mi/m-i, where the inverse-energy weighted mo- 
ment to_i is directly related to the static response [ p"2| . 
In the limit of a small harmonic perturbation Xr 2 , the 
induced density fluctuation takes the analytic form 



Sn = 



A / sin(7rr/i?o) 
27ru \ t/ttRq 



- 3 0{Rq - r) 



(16) 



and the inverse-energy weighted moment becomes to_i = 
-S(r 2 )/2\ = (6tt 3 - (2/5)7rf)m-i U -5?i 5 ai In terms 
of the plasma frequency (|l5|), we finally find lum = 
0.6181 ujp, a value slightly smaller than the one derived 
from the ratio 7713/7774. 

(b) G-regime — We turn now to the asymptotic G- 
rcgimc (self-bound "gravitational" regime) (T|, in which 
the s-wave interaction and the external trap potential 
are negligible, and "gravity" is balanced by the kinetic 
energy i/kin- In this case the monopole and quadrupole 
frequencies (|ll"l), (|l^ ) are given by wm — 0.3022 u p , ujq = 
0.5057 ujp, respectively, where 



lu p = 0.2351 mu 2 N 2 /h 3 



(17) 



is the plasma frequency (||) evaluated at the central den- 
sity in the same regime. The ratio between these frequen- 
cies is now u>m/wq = a/5/14 < 1. The chemical poten- 
tial \i = (i/kin + Hho + 2-ffgrav + 2H S ) / N in this regime 
has the value /1 = —0.6922 hw v . This shows that both 
the monopole and quadrupole oscillation frequencies are 
smaller than —fi/h. Consequently, these oscillations are 
stable against the ejection of atoms from the condensate. 

Also in the G-regime the evaluation of the static 
monopole response can be used to estimate the monopole 



frequency through the ratio Ti 



2 , ,2 



?77i/?77_i. The calcu- 



lation yields ujm — 0.2928 u> p . Similarly to the the TF-G 
regime the monopole frequencies evaluated through the 
ratios 7773/7771 and m\jm-\ are remarkably close, con- 
firming the validity of the method. 

CoUisionless hydrodynamics — In the Thomas-Fermi 
(TF-G) regime, exact results for the oscillation frequen- 
cies are obtainable by solving explicitly the equations of 
coUisionless hydrodynamics (|[), (||). These results allow 
us to improve the sum-rule estimate, which in general 
provides only an upper bound. The equations of mo- 
tion (|5|) and (^|) can be simplified in this regime, since 
the quantum pressure term, proportional to h 2 , can be 
safely ignored. After linearization they yield the equation 
of motion for the density fluctuations Sn 
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dt 2 



5n = V [c 2 s V (Sn + 6$ / g)] 



(18) 



where the "local" sound velocity c s (r) = y/gn(r)/m is 
calculated at the equilibrium density (0). The time de- 
pendence of Sn is chosen in the form of exp(— iuit), and 
(Eq) is transformed, with the help of (f|), into 



w 2 V 2 £ + V 



sin(27rf) - 
2irr 




= 



(19) 



where we have introduced the " gravitational" fluctuation 
variable £ = S<& / g and used the dimcnsionlcss coordinate 
f = r/Aj and frequency Cj — lu/u> p , with uj p given by 

To solve this equation, we set £ = ^x(f)Yi m i^ , 9 s ) ; where 
Yi m are the usual spherical harmonics. The boundary 
conditions can be derived from the tailoring of the inner 
and outer solutions on the condensate surface, i.e. at 
f = 7r jL4|. The conservation of the total number N of 
atoms implies that the radial derivative of the potential 
perturbation vanishes at r < n for the monopole mode. 
For the modes with I = 1,2, 3, ... the function instead 
decreases as r~~ Hence, the boundary conditions are 
(d4o/df)\ f=7r = for I = and (dft/df)^ = -(/ + 
l)f,/7r for i > 0. 

The numerical solutions of ( |l9|) for the monopole and 
quadrupole frequencies yield wm = 0.6168 u> p , ujq 



0.392 uj p with ui p as in (|15|). Hence, the sum-rule esti- 
mates given above, for these frequencies turn out to be 
very close (within 1%) to the exact numerical values. 
Furthermore, they satisfy the inequalities predicted by 
the general formalism of sum rules |Q . 

Conclusions — In this paper we have mainly dealt 
with the oscillations of a "gravitationally" self-bound 
bosonic cloud, in the framework of the (generalized) 
Gross-Pitaevskii approach. The Jeans wavelength (j^), 
which is the shortest wavelength to ensure stability, has 
been shown to be the diameter of the self-bound cloud 
in the Thomas-Fermi limit [Eq. @]. We have derived 
here simple, analytical expressions for the frequencies of 
the lowest monopole and quadrupole modes, which are 
the most relevant collective excitations of such a system, 
through a sum-rule approach. We have also numerically 
calculated the eigenfrequencies of the linearized general- 
ized Gross-Pitaevskii equation in the Thomas- Fermi limit 
of a self-bound BEC and demonstrated that the results 
are very close to the sum-rule values. 

The studied excitations substantially differ from those 
of harmonically trapped gases interacting via zero-range 
forces HfjJI]. For harmonically-trapped bosons, all the 
low-energy excitations (monopole, dipole, quadrupole 
etc.) are fixed by the oscillator frequency ujq of the confin- 
ing trap. By contrast, in the self-bound regimes, only the 
frequency of the lowest dipole mode, associated with the 
center of mass oscillations, is fixed by the external trap 
frequency ujq, whereas all the other low-energy excita- 
tions arc dominated by the effective "plasma" frequency 
w» (@). 



Since the repulsive short-range (s-wave) force can dras- 
tically change the size of the self-bound atomic cloud [[l] , 
the explicit values of u p and the corresponding monopole 
and quadrupole frequencies are determined by the central 
density and, consequently, by the regime considered [Eqs. 
(|l7| ) and (|l5|)1 . The dependence of the low-frequency col- 
lective modes on the plasma frequency resembles the os- 
cillations of compact stars, such as white dwarfs and neu- 
tron stars, where the frequency scale tu p ~ AirGp provides 
the universal relation between the oscillation timescale 
and the mass density p [Q. 
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